Entanglement of atoms via cold controlled collisions 
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We show that by using cold controlled collisions between 
two atoms one can achieve conditional dynamics in mov- 
ing trap potentials. We discuss implementing two qubit 
quantum-gates and efficient creation of highly entangled 
states of many atoms in optical lattices. 

PACS: 03.67.-a, 32.80.Pj, 03.67.Lx, 34.90.+q. 



The controlled manipulation of entangled states of N- 
particle systems is fundamental to the study of basic as- 
pects of quantum theory |l],|| , and provides the basis of 
applications such as quantum computing and quantum 
communications (^Q). Engineering entanglement in real 
physical systems requires precise control of the Hamilto- 
nian operations and a high degree of coherence. Achiev- 
ing these conditions is extremely demanding, and only 
a few systems, including trapped ions, cavity QED and 
NMR, have been identified as possible candidates to im- 
plement quantum logic in the laboratory ||. On the 
other hand, in atomic physics with neutral atoms recent 
advances in cooling and trapping have led to an exciting 
new generation of experiments with Bose condensates || , 
experiments with optical lattices Q , and atom optics and 
interferometry. The question therefore arises, to what ex- 
tent these new experimental possibilities and the under- 
lying physics can be adapted to provide a new perspective 
in the field of experimental quantum computing. 

In this Letter we propose coherent cold collisions as 
the basic mechanism to entangle neutral atoms. The 
picture of atomic collisions as coherent interactions has 
emerged during the last few years in the studies of Bose 
Einstein condensation (BEC) of ultracold gases. In a 
field theoretic language these interactions correspond to 
Hamiltonians which are quartic in the atomic field oper- 
ators, analogous to Kerr nonlinearities between photons 
in quantum optics. By storing ultracold atoms in arrays 
of microscopic potentials provided, for example, by opti- 
cal lattices these collisional interactions can be controlled 
via laser parameters. Furthermore, these nonlinear atom- 
atom interactions can be large even for interactions 
between individual pairs of atoms, thus providing the 
necessary ingredients to implement quantum logic. 

Let us consider a situation where two atoms \a) and 
\b) are trapped in the ground states ip ' of two poten- 
tial wells V a ' b . Initially, at time t = — r, these wells 
are centered at positions x a and x b , sufficiently far apart 
(distance d = Xb — x a ) so that the particles do not in- 



teract. The positions of the potentials are moved along 
trajectories x a (t) and x b (t) so that the wavepackets of 
the atoms overlap for certain time, until finally they are 
restored to the initial position at the final time t = r. 
This situation is described by the Hamiltonian 
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and p a ' b are position and momentum opera- 
(x a ' b — x a ' b (t)) describe the displaced trap po- 
tentials and u ab is the atom-atom interaction term. 

Ideally, we would like to implement the transformation 
from t — — t to t = t 

%{x a -x a )^(x b -x b ) -> e l(t >^{x a -x a )i) b {x b -x b ), (2) 

where each atom remains in the ground state of its trap- 
ping potential and and preserves its internal state. The 
phase <fi will contain a contribution from the interac- 
tion (collision). Transformation (^) can be realized in 
the adiabatic limit, || whereby we move the potentials 
so that the atoms remain in the ground state. In the 
absence of interactions (u ab = 0) adiabaticity requires 

\x a ' <C v OBC Vt, where v OS c ~ oquj is the rms veloc- 
ity of the atoms in the vibrational ground state, ao is 
the size of the ground state of the trap potential, and ui 
is the excitation frequency. The phase </> can be easily 
calculated in the limit \x a ' (t)\ <C v osc /t. In this case, 
6 — 6 a + 6 b where 
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are the kinetic phases. In the presence of interactions 
(w ab 7^ 0), we define the time-dependent energy shift 
due to the interaction as 
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where a s is the s-wave scattering length. We assume 
that: (i) |Ai5(i)| <C hui so that no sloshing motion is ex- 
cited; (ii) |a; a ' b (£)| <ti v osc (adiabatic condition); (iii) u osc 
is sufficiently small for the zero energy s-wave scattering 
approximation to be valid ||. In that case, (||) still holds 
with 4> = 4> a + <f) b + ab , where 
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is the collisional phase. 

In the case of quasi-harmonic potentials (as is realized 
with optical potentials, see below) (||) still holds, even in 
the non-adiabatic regime, i.e. at higher velocities. For 
harmonic traps one can solve exactly the evolution for 



, ab 



= 0, and identify the condition for adiabaticity: 
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This is consistent with condition (ii) . Actually, for (||) to 
hold the inequality (JsJ) only has to be fulfilled for t = t 
(and not for all times). This means that the particle 
need not be in the ground state of the moving poten- 
tial at all times, but only at the final time. The phase 
(j) = <p a + 4> b , as well as the wave functions ip a ' b (x a ' b , t), 
can be calculated exactly ||. For u ab ^= if conditions 
(i) and (iii) are satisfied, then (||) is still valid. There 
is an additional phaseshift (with respect to u ab = 0) 
</> ab which is given by Eqs. with the replacement 

ii%' b (x a > b - x a - b (t)) -► iP a < b (x a ' b ,t). It is also straightfor- 
ward to generalize these results to the case in which the 
trap frequency changes with time JTo| . 

So far, we have shown that one can use cold collisions 
as a coherent mechanism to induce phase shifts in two- 
atom interactions in a controlled way. Our goal is now 
to use these interactions to implement conditional dy- 
namics. We consider two atoms 1 and 2, each of them 
with two internal levels |a)i ; 2 and |6)i,2- We will use the 
superscripts /3 = a,b and the subscripts j = 1,2 to la- 
bel the internal levels and atoms, respectively. Atoms in 
the internal state \/3)j experience a potential which 
is initially (t = — r) centered at position Xj. We assume 
that we can move the centers of the potentials as follows 
(Fig. 1): x?(t) = Xj+5xP(t). The trajectories 5x (t) are 
chosen in such a way that 8xP(-t) = 5xP{t) = and the 
first atom collides with the second one only if they are in 
states | a) and \b), respectively (\x\(t) — x 2 {t)\ ^> ao Vi). 
This choice is motivated by the physical implementation 
considered below. The fact that Xj does not depend on 
the internal atomic state allows one to easily change this 
internal state at times t = ±r by applying laser pulses. 
If the conditions stated above are fulfilled, depending on 
the initial internal atomic states we have the following 
table for this process: 

|a) 1 |a) 2 -> e ' i2 * a |a) 1 |a) 2 , 

\a)i\b) 2 ^e^'+* b+ *^\ah\b) 2 , 

\b) 1 \a) 2 ^e^ a +^\b) 1 \a) 2 , 

|fo)i|fo) 2 -e i20 >)i|6) 2 , 
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where the motional states remain unchanged. The ki- 
netic phases <jr and the collisional phase </> ab can be cal- 
culated as stated above. We emphasize that the <fiP are 
(trivial) one particle phases that, if known, can always 
be incorporated in the definition of the states \a) and \b). 




FIG. 1. Configurations at times ±r (a) and at t (b). The 
solid (dashed) curves show the potentials for particles in the 
internal state |o) (|6)), respectively. Center positions (t) 
and displacements Sx^(i) as defined in the text. 

In the language of Quantum Information, transforma- 
tion (R) corresponds to a fundamental two-qubit gate 
[Q. In order to illustrate to what extent the mechanism 
presented above is able to perform this ideal gate, we 
have carried out a numerical study in 3 dimensions. We 
have integrated the time-dependent Schrodinger equa- 
tion with the Hamiltonian (Q). We have taken harmonic 
potentials with various time dependent displacements 
5x@(t) and frequencies. Their form as well as the param- 
eter range are motivated by the specific implementations 
outlined below. The figure of merit that we have used is 
the minimum fidelity, which is defined as 



F = min<#r ext (U\ip){ip\ <8> Pcxt U^ 
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Here \ip) is an arbitrary internal state of both atoms, 
is the state resulting from \ip) using the mapping <^). 
The trace is taken over motional states, U is the evolu- 
tion operator for the internal states coupled to the exter- 
nal motion (including the collision), and p c ^ the density 
operator corresponding to both atoms in the motional 
ground state at t = — r. In the ideal case the fidelity 
will be one. We have taken the potential u ab as propor- 
tional to a delta function and used a truncated moving 
harmonic oscillator basis (10 states for each degree of 
freedom). The first illustration corresponds to a fixed 
trap frequency io a = ui b and displacements 5x b (t) = 
and 8x a (t) as specified in Fig. 2a. In Fig. 2b we present 
a contour plot of F as a function of the parameters char- 
acterizing the displacement 5x a (t). The fidelity is very 
close to one for a surprisingly wide range of parameters, 
even well into the non-adiabatic regime. We have also 
studied numerically the effect of time-varying trap fre- 
quencies (see Fig. 3a) and finite temperatures, so that 
/9 ex t describes a thermal distribution of temperature T. 
For temperatures kT < 0.2hui (where u) is the initial trap 
frequency) the fidelity remains very close to one (of the 
order of 0.997), and this result is still robust with respect 
to changes in the parameters. We have also included the 
possibility of loss due to collisions by adding an imag- 
inary part to the scattering length. In that case, the 
fidelity is reduced (see Fig. 3b). 
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FIG. 2. a) 

8x a (t)/d = (l+exp(-(r l /r r ) 2 ))/(l+e X p((t 2 -r/)/r r 2 )) 
and 5x b (t) = (see text), with r r — 30/u; and Ti = 20/oj. 
The shaded region indicates where the particles interact, b) 
Fidelity F against rise time r r and interaction time n for 
87 Rb with a s — 5.1nm, uj = 2tt x 100kHz and d — 10 do. 

So far we have assumed that there is one atom per 
potential well. In practice, the particle number might 
not be controlled. Nevertheless, one can easily generalize 
the above results to this case by using a second quantized 
picture. For example, in the adiabatic regime we denote 
by a,i and bi the annihilation operators for a particle in 
the ground state of the potential centered at the position 
i, and corresponding to the internal levels \a) and \b), 
respectively. The effective Hamiltonian in this regime is 

H = Y^ [^ a (t)4 a i + u b (t)bjb t + u aa {t)a\a\aiai+ 

i 

u hh (t)b\b\bA] +Y.u^{t)a\ ai b%, (9) 

where the w's and u's depend on the specific way the 
potentials are moved. This Hamiltonian corresponds to 
a Quantum-Non-Demolition situation fPH , whereby the 
particle number can be measured non-destructively. 




FIG. 3. a) Upper plot: Displacements 8x a (t)/d (solid 
line) and 1 + 5x b (t)/d (dashed line). Lower plot: Trap 
frequencies uj a {t)/ui (solid line) and uj (t)/uJ (dashed line), 
b) Fidelity F against temperature kT/Tiuj for 87 Rb with 
a 3 = 5.1nm (solid line), a s — (1 — O.Oli) x 5.1nm (dashed 
line) and a 3 = (1 — 0.05i) x 5.1nm (dash-dotted line). Here 
uj = 2tt x 100kHz and d = 390nm. 



A physical implementation of this scenario requires an 
interaction which produces internal state dependent con- 
servative trap potentials and the possibility of moving 
these potentials independently. Furthermore, the choice 
of the internal atomic states |o) and \b) has to be such 
that they are collisionally stable (i.e. the internal states 
do not change after the collision). These requirements 



can be satisfied in an optical lattice. We consider a spe- 
cific but relevant example of alkali atoms with a nuclear 
spin equal to 3/2 ( 87 Rb, 23 Na) trapped by standing waves 
in 3 dimensions. The internal states of interest are hyper- 
fine levels corresponding to the ground state Si/ 2 . Along 
the z axis, the standing waves are in the linZlin configu- 
ration (two linearly polarized counter-propagating trav- 
eling waves with the electric fields forming an angle 29 
p"2|). The electric field is a superposition of right and 
left circular polarized standing waves (cr ) which can be 
shifted with respect to each other by changing 9, 

E + (z,t) = E e- wt [e + sin(kz + 9) + e_sm(kz-9)}, (10) 

where e± denote unit right and left circular polarization 
vectors, k — v/c is the laser wavevector and Eq the ampli- 
tude. The lasers are tuned between the P\/ 2 and P 3 / 2 lev- 
els so that the dynamical polarizabilities of the two fine 
structure S1/2 states corresponding to m s = ±1/2 due to 
the laser polarization vanish, whereas the ones due to 
a are identical (= a). The optical potentials for these 
two states are V ma= ±i/ 2 (z, 9) = a|i?o| 2 sin 2 (kz ± 9). We 
choose for the states \a) and |6) the hyperfme structure 
states I a) = \F = l,m/ = 1) and \b) = \F = 2,m f = 2). 
Due to angular momentum conservation, these states are 
stable under collisions (for the dominant central elec- 
tronic interaction The potentials "seen" by the 
atoms in these internal states are 

V a (z, 9) = [V ms=1/2 (z, 9) + 3V ms= ^ 1/2 ( Zl 9)] /4 (11a) 



V\z,9) = V ms=1/2 (z,9). 



(lib) 



If one stores atoms in these potentials and they are deep 
enough, there is no tunneling to neighboring wells and 
we can approximate them by harmonic potentials. By 
varying the angle 9 from 7r/2 to 0, the potentials V b 
and V a move in opposite directions until they completely 
overlap. Then, going back to 9 = tt/2 the potentials re- 
turn to their original positions. The shape of the po- 
tential V a changes as it moves. By choosing 9(t) = 
tt (1 - (1 + exp(-(r l /r r ) 2 )) / (l + exp((* 2 - r 2 )/r 2 ))) /2 
with T r = 30/w and r, = 20 /to, the frequencies and 
displacements of the harmonic potentials approximating 
( pT| ) are exactly those plotted in Fig. 3a. Therefore, that 
figure shows that under this realistic situation one can 
obtain very high fidelities. 

The scheme presented here can be used in several in- 
teresting experiments: (a) One can use this method to 
measure the phase shift <fi ah and therefore determine the 
scattering length corresponding to a— b collisions. For 
that, one can use ideas borrowed from Ramsey spec- 
troscopy, namely: (i) prepare all atoms in the superpo- 
sition (| a) + |b))/V2 by applying a tt/2 laser pulse; (ii) 
shift their potentials as described above; (iii) apply an- 
other tt/2 laser pulse; (iv) detect the population of the 
internal states by fluorescence. It can easily be shown 
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that such populations depend in a simple way on the 
phase shift ab . One can even determine the sign of 
the scattering length by applying laser pulses with pulse 
area different from tt/2 [l^] . (b) In a similar way, one 
can also measure the spatial correlation function by ap- 
plying the second laser pulse (iii) and population detec- 
tion (iv) without moving the potential back to the origin. 
This would be a way of discriminating between Mott and 
superfluid phases for particles in an optical lattice [Q. 
(c) Apart from that, if one is able to address individual 
wells with a laser, one can also perform certain experi- 
ments which are interesting both from the quantum in- 
formation and fundamental point of view. For example, 
one could create an entangled EPR pair of two particles 
(|a)i|6) 2 -|6)i|a) 2 )/V2 |]. This could be done by having 
two atoms in neighboring wells and: (i) prepare each of 
them in the superposition (|a) + \b))/\/2; (ii) shift their 
potentials back and forth so that the phase shift cf) ah = n: 
(iii) apply a ir/2 laser pulse to the second atom. In this 
case one could test Bell inequalities and perform other 
fundamental experiments, (d) With more than two par- 
ticles one could create higher entangled states with sim- 
ple lattice operations. For example, if one has N parti- 
cles in N potential wells, one could create GHZ states of 
the form (|o)i|o) 2 . . . | a) jv - | &) 1 1 &) 2 . . . \b) N )/y/2 @. This 
could be easily done, for example, if one could use a dif- 
ferent internal state |c)i in the first atom instead of 
jl5||. The procedure would be as follows: (i) prepare the 
first atom in the state (|a)i + |c)i)/\/2, and the others 
in the state (|o) + \b))/y/2; (ii) Move the potential corre- 
sponding to the internal level |c) back and forth so that 
if the first atom is in that state it interacts with all the 
other atoms for a time such that in each "collision" the 
phase shift difference between the collision a-c and b-c 
is 7r; (iii) apply a n/2 pulse (in transition a—b) to all the 
atoms except the first one, which is transferred from \ c) to 
1 6). We emphasize that the optical lattice configuration 
allows the preparation of these states in a single sweep of 
the lattice, (e) Finally, it is clear that one could use op- 
tical lattices in the context of quantum information since 
the above procedure provides a fundamental two bit gate 
(0) which, combined with single particle rotations allows 
to perform any quantum computation between an arbi- 
trary number of two-level systems. In particular, the 
optical lattice setup would be very well suited to imple- 
ment fault tolerant quantum computations due to the 
possibility of doing gate operations m in parallel (details 
will be given in [H). 

So far we have neglected some processes that may lead 
to decoherence, and therefore limit the performance of 
our scheme. This includes spontaneous emission of atoms 
in the off-resonant optical lattice potentials, and inelastic 
collisions to wrong final atomic states. The spontaneous 
emission lifetime of a single atom in the lattice is from 
seconds to many minutes, depending on the laser detun- 
ing |6j . The problem of collisional loss is closely related to 



the loss mechanism in Bose condensates in magnetic and 
optical traps, a problem studied extensively in recent ex- 
periments. By choosing proper internal hyperfine states 
one can maximize the lifetime due to inelastic collisions 
to be of the order of at least several seconds jl6| . 

Finally, by filling the lattice from a Bose condensate, 
and using the ideas related to Mott transitions in optical 
lattices H it is possible to achieve uniform lattice occupa- 
tion ("optical crystals") or even specific atomic patterns, 
as well as the low temperatures necessary for performing 
the experiments proposed in this letter. 

This research was supported by the Austrian Science 
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Note added: After this work was completed we became 
aware of G. K. Brennen et al, quant- ph/9806021 , where 
dipole-dipole interactions are proposed as a mechanism 
to implement a two bit quantum gate. 
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